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1.  Introduction:  The  Nonlinear  Transmission  Line 


Transmission  lines  are  ubiquitous  in  almost  all  high-frequency  RF  circuits,  where  they  replace 
the  standard  wire  connections  used  for  assembling  circuits  at  lower  frequencies.  The  many  forms 
they  take,  involving  complex  geometries  and  exotic  materials,  are  conditioned  by  the  task  of 
delivering  RF  energy  from  one  circuit  location  to  another.  A  less  common  use  of  transmission 
lines  is  as  actual  RF  devices  -  switches,  phase  shifters,  filters,  etc.  -  which  exploit  their  ability  to 
perform  various  impedance  functions  and  versatility  with  regard  to  fabrication.  For  example, 
varying  lengths  of  the  same  transmission  line  can  be  used  as  capacitors,  inductors,  or  delay  lines 
for  matching  circuit  subunits  such  as  amplifiers  (1).  When  used  in  this  way,  the  transmission  line 
is  usually  a  linear  circuit  component.  A  considerably  less  common  use  of  nonlinear  transmission 
lines  (NLTL)  is  for  large-signal  device  design  -  phase  shifters  come  to  mind  (2)  -  although  some 
early  work  on  distributed  amplifiers  could  qualify  as  such  (5). 

The  structure  of  a  transmission  line  can  be  either  continuous  or  discrete.  Whereas  continuous 
transmission  lines  are  typically  waveguide-based  (5a),  lines  with  discrete  structure  usually 
consist  of  repeated  identical  subcircuits  made  up  of  lumped  elements.  A  useful  definition  of  a 
NLTL  is  an  ordinary  transmission  line  whose  electrical  characteristics  are  functions  of  local 
voltages  and  currents.  Thus,  the  prototypical  discrete  transmission  line  shown  in  figure  1  is  made 
up  of  lumped  elements  -  resistors,  capacitors,  inductors,  etc.  -  any  one  of  which  can  be  voltage 
dependent.  In  the  1980’s,  NLTL  were  fabricated  by  inserting  shunt  diodes  periodically  into  a 
transmission  line  to  create  nonlinearity  in  its  overall  electrical  response  (3b).  Although  these 
lines  had  very  useful  properties,  they  were  never  capable  of  much  power  handling  due  to  the 
semiconductor  nature  of  the  diodes  used.  At  the  end  of  the  1980’s,  the  Ikezi  group  at  General 
Atomics  (4)  studied  layered  structures  consisting  of  parallel-plate  waveguides  in  which  layers  of 
ferroelectric  material  (specifically,  barium  strontium  titanate,  or  BST)  alternated  with  layers  of 
simple  dielectric  (see  figure  2);  these  structures  now  represent  the  canonical  version  of  high- 
voltage  NLTLs,  and  are  often  encountered  in  the  NLTL  literature  to  this  day.  Mark  Litz  and  I 
were  first  made  aware  of  the  work  of  Ikezi’s  group  in  2007,  as  a  possible  way  to  address  the  need 
of  ARDEC  for  a  HPM  source,  and  it  has  been  the  basis  for  our  effort  at  ARL  ever  since. 

In  this  report  I  will  describe  my  efforts  to  model  the  electromagnetic  response  of  this  type  of 
transmission  line  numerically,  as  requested  by  ARDEC  in  the  summer  of  2009. 
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Figure  2.  NLTL  structure  used  for  fabrication  at  ARL. 


2.  Selection  of  Capacitance- Voltage  Dependence  for  NLTL  Circuit 


Much  of  the  NLTL  design  work  published  in  the  literature  is  numerieal  in  nature.  The  need  to 
use  numerical  modeling  to  design  a  useful  NLTL  is  regrettable  in  many  ways,  the  most  important 
of  which  is  that  easy-to-use  design  rules  are  seldom  derivable  from  simulations.  After  much 
labor,  Ikezi  et  ah,  at  General  Atomics  (5)  published  an  example  of  one  sueh  design  rule  for 
finding  the  frequency  at  which  the  power  from  a  NLTL  was  a  maximum.  At  ARL  we  used  this 
rule  extensively  in  our  previous  effort.  The  design  rule  was  formulated  by  performing  a  series  of 
numerical  simulations  of  their  NLTL  eireuit  model  for  various  line  parameters,  and  extraeting 
the  peak  power  frequency  by  Fourier- transforming  the  computed  time  dependence  of  the  power 
out  of  the  line  and  looking  for  its  peak  value.  Their  results  were  published  in  the  form  of  a  curve, 
whieh  they  fitted  to  the  empirical  function 


—  =  1.5  logio 


+  1.25 
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where  co  is  the  peak-power  frequency  of  the  output  soliton  pulse  train,  co^  is  the  width  of  the 
lowest  NLTL  pass  band  when  treated  as  a  linear  (small-signal)  transmission  line,  E  is  the  electric 
field  measured  at  the  input  of  the  NLTL,  and  Eq  is  a  scaling  field  of  order  10  kV.  Because  the 
relationship  implied  by  this  formula  between  the  small-signal  and  large-signal  quantities  is 
obscure,  its  usefulness  for  evaluating  other  systems  is  problematic.  Because  we  were  unable  to 
make  any  accurate  power  measurements,  we  were  unable  to  check  how  good  the  equation  was; 
nevertheless,  we  were  able  to  use  the  dependence  on  co^  to  design  the  NLTL  at  a  prespecified 
frequency,  using  an  additional  equation  taken  from  (6). 

Because  the  nonlinear  response  of  the  NLTL  is  based  on  the  voltage  dependence  of  BST 
capacitors,  it  is  vital  that  these  components  be  modeled  as  precisely  as  possible.  At  this  time  we 
have  two  literature -based  modeling  functions,  one  the  original  choice  of  Ikezi  in  1988  (5)  and  the 
other  that  of  Martin  Brown  in  his  1997  thesis  (7).  In  this  section  I  will  explain  why  the  former 
was  preferred  over  the  latter  as  a  basis  for  this  modeling  work,  and  propose  a  third  alternative 
with  features  of  both. 

2,1  Ikezi’s  Model 

The  model  used  by  the  Ikezi  group  for  the  voltage  dependence  of  BST  capacitors  is  based  on  a 
fundamental  relationship  between  the  electric  displacement  in  BST  and  an  external  electric  field, 
which  they  write  as  follows: 

^  D 

E  - — T — - 

The  quantity  D  gives  the  charge  per  unit  area  induced  on  the  metal  walls  of  the  parallel-plate 
transmission  line,  and  hence,  can  be  multiplied  by  the  area  of  contact  A  between  the  metal  and  a 
given  ferroelectric  layer  to  extract  the  charge  per  ferroelectric  layer  induced  on  the  walls: 

D  =  ^^Q  =  AD 
A 

while  the  electric  field  E  equals  the  voltage  F  between  the  walls  divided  by  the  height  of  the 
guide: 

V 

E  =  —  ^V  =  Eh 
h 

Then  together  these  quantities  let  us  define  a  function  V{Q)  which  can  be  used  to  specify  the 
capacitance  of  the  line: 
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Note  that  in  order  to  eompute  the  eapaeitanee,  we  must  invert  the  funetion  V(Q)  numerieally. 

The  form  of  this  expression,  whieh  is  a  eonsequenee  of  the  erystal  physies  of  BST  (5),  predicts  a 
capacitance  versus  voltage  with  the  following  properties: 

1 .  At  small  voltages  (small  D)  it  predicts  a  constant  capacitance  of  the  usual  form,  i.e., 

r  _  ^^0^ 

h  ’ 

2.  The  capacitance  is  independent  of  polarity,  and 

3.  For  large  voltages  the  capacitance  decreases  algebraically  as  V 

It  is  worth  noting  that  the  expression  used  by  Ikezi  is  clearly  a  truncation  of  an  infinite  series, 
and  indeed,  it  is  common  practice  with  in  most  of  the  ferroelectrics  literature  to  write 

E  =  qD  +  bD  +  cE) 

-4/5 

For  our  purposes,  it  is  significant  that  this  implies  a  falloff  in  the  capacitance  as  V  with 
voltage.  The  large- F behavior  is  thus  ambiguous,  and  can  only  be  resolved  by  fitting.  Note  that 
the  procedure  for  fitting  experimental  data  to  polynomials  is  especially  simple,  a  further 
advantage  to  using  Ikezi ’s  expressions. 

2,2  Brown’s  Model 

In  contrast  to  the  Ikezi  physics-based  model,  the  model  used  by  Brown  seems  to  be  entirely  ad 
hoc,  and  violates  both  the  material  constraints  and  the  standard  laws  of  electrostatics,  casting  his 
entire  simulation  effort  into  doubt: 

(1)  He  assumes  that  the  relationship  between  charge  and  voltage  is 

Q  =  C{V)V 
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whereas  in  fact  the  correct  relationship  is 


e(F)=jc(u)^/u 


(2)  He  postulates  a  diodelike  functional  dependence  for  C(V)  that  depends  on  polarity: 

C(H)  =  Cora  +  (l-a)e“''^^" 

In  addition,  this  expression  predicts  a  strong  increase  for  negative  F,  which  is  unphysical. 

It  is  interesting  to  determine  the  difference  between  Brown’s  expression  and  the  correct 
expression  for  Q(V):  his  simple  multiplication  gives 


Q(F)  =  C(F)F  =  CoFf  a  +  (l-a)e 


-V  I  Vr 


while  the  correct  expression  gives 


V 

g(H)  =  Co|[a  +  (l-a)e 


-V  IV r 


dV  =  C<, 


a 


V  +  {l-a)Vc{l-e 


^-V  I  Vr 


The  plot  shown  in  figure  3  illustrates  the  difference: 


Figure  3.  The  function  Q(V)  evaluated  from  Brown’s  fitting  expression:  red  -  multiplication, 
blue  -  integration. 
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Clearly  the  two  functions  are  almost  identical  for  small  voltages,  and  differ  only  by  an  offset  for 
large  voltages;  however,  the  severe  change  in  slope  at  intermediate  voltages  will  greatly  distort 
the  waveform  of  an  ac  signal  as  it  passes  through  the  intermediate  voltage  range. 

These  deficiencies  persuade  us  to  adopt  the  original  Ikezi  relations  in  our  simulations. 


3.  Nonlinear  Circuit  Software  Development  For  Transmission  Line 


Because  of  its  flexibility  and  proven  track  record  with  regard  to  computational  accuracy  and 
speed,  the  mathematical  package  Mathematica™  has  been  chosen  for  software  support  of  the 
NLTL  design.  In  order  to  reduce  the  learning-curve  burden  for  future  users,  code  written  for  this 
package  will  be  translated  into  MathCad™  in  2010  if  funding  permits;  meanwhile,  the  basic 
analysis  of  the  NLTL  performance  benefits  greatly  from  the  feasibility  of  writing  transparent 
software  that  is  easily  changed  and  physically  motivated,  and  is  tailored  to  the  specifics  of  the 
NLTL  problem  rather  than  relying  on  more  broadly  applicable  circuit  software  such  as  PSPICE 
which  cannot  be  “tuned”  internally,  i.e.,  lack  of  access  to  source  code. 

The  approach  to  circuit  modeling  adopted  here  takes  advantage  of  the  simplicity  of  the  NLTL 
circuit  model  which  in  turn  is  a  consequence  of  the  underlying  parallel-plate  waveguide 
geometry  of  the  latter.  A  more  complex  propagation  structure  would  require  much  more 
powerful  mathematical  methods  to  treat  correctly.  The  modeling  will  be  described  in  the 
following  sections. 

3,1  Brown  and  Ikezi  Circuits 

Both  Brown  and  Ikezi  adopt  a  simple  series-parallel  circuit  model  for  the  NLTL  (see  figure  4). 


A  justification  of  this  circuit  structure  will  be  given  in  a  future  report.  The  original  circuit 
analyzed  at  General  Atomics  had  no  resistances;  the  various  resistors  shown  here  were  added  by 
Brown.  The  circuit  is  an  oversimplification  of  the  real  NLTL;  roughly  speaking,  the  inductors 
represent  the  alumina  layers,  while  the  shunt  lines  represent  the  BST  layers.  This 
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oversimplification  will  eventually  be  removed  in  our  modeling  effort.  The  capacitors  C  are  taken 
to  be  nonlinear  elements. 


3.2  State  Variables 

Before  the  advent  of  circuit  simulators,  it  was  well  understood  that  nonlinear  circuits  must  be 
approached  as  individual  problems.  This  is  in  contrast  to  linear  circuits,  which  allow  a  universal 
approach  to  computations  based  on  straightforward  matrix  computations.  Textbooks  on  circuit 
theory  routinely  introduce  the  idea  of  the  “state”  of  a  circuit,  by  which  they  mean  a  list  of  all  the 
currents  flowing  through  the  various  circuit  components  and  all  the  voltage  drops  across  them. 
When  these  quantities  depend  on  time,  as  in  an  ac  circuit,  their  values  are  obtained  from 
differential  equations.  For  the  special  case  of  linear  circuits  driven  by  sinusoidal  sources,  these 
differential  equations  reduce  to  algebraic  equations  that  can  easily  be  solved  by  matrix  methods. 
Transient  excitation,  e.g.,  by  pulses,  require  the  original  differential  equations  and  numerical 
methods,  as  do  circuits  with  nonlinear  elements. 

It  is  important  for  numerical  reasons  that  these  differential  equations  all  be  first-order  in  time, 
i.e.,  if  we  describe  the  circuit  state  using  a  vector  of  currents  and  voltages 
f  =  (  Fj ,  ,  F2 ,  /2 ,  F3 ,  /j  ,•••),  the  numerical  problem  must  be  of  the  form 


dx 

dt 


=  F(x) 


(1) 


where  F’(x)  is  a  vector  function  of  the  vector  variable  x  .  Finding  this  function  can  be  very 
difficult  for  a  complicated  circuit  topology,  since  the  circuit  equations  usually  take  the  form 


F 


X, 


dx 

dt 


=  0 


(2) 


Numerical  simulators  like  PSPICE  routinely  use  a  numerical  engine  (a  “root  solver”)  that 
converts  equation  2  into  equation  1,  which  greatly  increases  execution  time  since  it  must  be  done 
at  each  time  step. 

It  is  usually  found  that  not  all  the  currents  and  voltages  are  needed  to  describe  the  circuit 
dynamics.  Thus,  for  a  simple  ac  circuit  like  the  one  in  figure  5;  only  two  differential  equations 
are  needed,  the  currents  through  the  two  inductors.  These  currents  are  said  to  be  the  state 
variables  of  the  circuit.  In  general,  state  variables  are  defined  as  the  smallest  set  of  electrical 
quantities  that  can  completely  specify  the  circuit  performance. 
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For  more  complex  circuits  there  may  be  multiple  choices  for  these  state  variables.  While  currents 
and  voltages  are  appropriate  for  resistive  elements,  the  natural  variables  for  capacitors  and 
inductors  are  actually  charges  and  magnetic  fluxes,  and  these  quantities  are  often  handy 
alternatives  to  currents  and  voltages,  especially  when  the  circuit  does  not  immediately  admit  a 
description  using  first-order  differential  equations. 

A  trivial  example  is  a  series  LRC  circuit  (see  figure  6). 


This  circuit  is  usually  treated  by  writing  a  single  equation 


Vo  =  Vl+V^+Vc 


(3) 
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where  the  voltage  drops  are  given  by 


V,=L- 


dl 

dt 


Vr=R1 


This  equation  ean  be  solved  easily;  however,  from  a  numerieal  standpoint  it  has  the  problem  that 
it  has  an  integral  in  it,  whieh  means  it  is  not  a  differential  equation.  Taking  the  derivative  if 
equation  3  eorreets  this: 


JTq  _  ^  d^l 
dt  dt^ 


udi  1  , 
+  R  —  +— / 
dt  C 


Although  this  is  easily  solved  analytieally,  it  is  not  first  order  and  henee  is  problematie  for 
numerieal  solution.  A  new  ehoiee  of  state  variables  fixes  this:  sinee  the  eurrent  is  related  to  the 
eharge  on  the  eapaeitor  by 


/ 


dt 


we  ean  write  the  following  system  of  equations: 


dQ^I 

dt 

dt  L  ^  L  LC 


We  now  have  a  state  veetor  x  = 


^0^ 


to  deseribe  the  cireuit,  and  a  veetor  funetion 
/ 


f 


F{x)  = 


-Vq--i — ^—Q 

\L  ^  L  LC  J 


whieh  a  simulator  ean  aeeept. 


It  is  fortunate  that  this  ehoiee  of  state  variables  also  works  for  the  NLTL.  The  appendix  eontains 
the  derivation  of  the  set  of  network  equations  for  the  Brown-Ikezi  eireuit  for  an  NLTL  with  an 
arbitrary  number  of  stages.  This  set  of  equations  has  been  eoded  into  Mathematica,  and  seems  to 
work  well  -  eireuits  with  over  a  thousand  stages  have  been  simulated  without  numerieal 
problems. 


3,3  The  Method  of  Partial  Charges 

In  writing  the  network  equations  for  this  eireuit,  it  is  useful  to  generalize  the  state  variables 
diseussed  above.  As  diseussed  in  the  appendix,  a  natural  way  to  treat  the  NLTL  eireuit  is  by  loop 
eurrents,  as  shown  in  the  eireuit  pieture  at  the  beginning  of  this  report.  Now,  the  diseussion  in  the 
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previous  section  suggests  that  we  introduce  a  “charge”  qn  associated  with  each  loop  current, 
which  we  define  as 

.  _dq„ 

"  dt 

This  identifies  a  set  of  state  variables  ,  ^2  ’  ' ' '  {h  ■>  h  ■>  h }  the 

circuit,  which  a  differential-equation  solver  can  handle  easily.  Once  we  have  the  solution  to  the 
network  equations,  we  can  get  the  node  voltages  from  the  obvious  relation 

Vn  =tc«  -4+1) 

where  Vcn  is  the  voltage  drop  across  the  nonlinear  capacitor.  As  discussed  in  Report  #1,  we  write 
this  voltage  in  the  form 

Va,=C-'Q„+aQ„^ 

where  the  capacitor  charge  Qn  -  qn  -  dn+i  mimics  the  behavior  of  the  shunt  current.  This  is  a  true 
charge,  arising  from  the  flow  of  current  into  and  out  of  the  shunt  branch. 

It  is  remarkable  that  when  this  decomposition  of  Qn  into  “partial  charges”  qn  is  made, 
preliminary  simulation  results  reveal  that  the  qn  are  unphysical  even  though  their  difference  is 
not.  When  a  NLTL  is  pulsed  by  a  finite  voltage  signal,  after  a  long  time  all  the  variables  qn  are 
nonzero  and  identical  in  value.  This  means  that  the  currents  and  true  charges  on  the  capacitors 
must  vanish,  as  they  should  in  a  dissipative  network  (i.e.,  one  with  resistors).  This  rather 
mysterious  result  has  no  impact  on  the  modeling  results,  but  it  is  an  interesting  consequence  of 
the  freedom  to  assign  state  variables  to  a  nonlinear  circuit. 


4.  Experimental  Fitting  of  Resistive  Losses  in  the  NLTL  Based  on  Loss 
Tangent  Data 


Accurate  determination  of  losses  in  ferroelectric  materials  such  as  BST  is  extremely  difficult, 
especially  over  a  wide  frequency  range.  Because  such  losses  have  a  major  impact  on  the 
performance  of  our  NLTL,  particular  care  was  taken  to  identify  a  reasonable  number  for  the 
single  parameter  we  use  to  characterize  these  losses  -  the  resistance  r  shown  in  the  circuit  (see 
figure  7). 
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Neither  Brown  nor  Ikezi  provide  any  information  about  this  number,  although  Brown  suggests 
that  he  used  loss-tangent  data  to  obtain  a  value  for  r  (note  that  he  uses  the  eireuit  shown  above, 
just  as  we  do).  Beeause  loss-tangent  data  does  not  direetly  translate  into  r,  a  proeedure  was 
developed  for  eonneeting  these  numbers,  whieh  will  be  deseribed.  All  the  simulations  performed 
will  be  based  on  numbers  obtained  from  this  proeedure. 


Frequeney  dependent  measurements  of  BST  parameters  usually  are  based  on  eapaeitanee 
measurements  made  on  fabrieated  eapaeitor  struetures.  This  proeedure  is  typified  by  the  work  of 
Rimai  et  ah,  (P),  who  used  eapaeitors  to  measure  dieleetrie  eonstant  and  loss  tangent  as  a 
funetion  of  frequeney.  They  fitted  their  data  to  the  following  funetional  form; 


8  =  Sr 


|l  +  Xo  “■ 


ycoa-co 


cOj  +  ycoG-co" 


-[Xo-Xoo] 


where  the  parameters  a>t,  Xq,  Xqq  a  are  used  to  fit  the  data.  The  values  of  these  parameters 
they  found  for  their  partieular  BST  samples  were: 

CO,  =1020  GHz 
Xo=1900 

X»  =6 

a  =  2.5(0, 


Now,  the  loss  tangent  is  defined  in  most  of  the  ferroeleetrie  literature  as  follows: 


tan5(co) 


Im8((o) 

Res(co) 
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so  its  values  are  obtained  by  finding 


Re8(co)  =  8o  Re 


co^  +  ycoa-co 


=  SoRe<! 


(1  +  Xo )  -  (1  +  Xoo )  +  +  Xoo ) 


-co^  -j(0O 


2  2 
CO^  -CO 


-i2 


2  2 
+  C0  G 


<(l  +  Xo)-®Ml  +  Xoo)  <-co^  +coV(l  +  xJ 


2  2 
CO^  -CO 


n2 


2  2 
+  C0  G 


and  likewise 


Im8(co)  =  8o  Im 


®c^(l  +  Xo)  +  (7®f^-®^)(l  +  Xoo) 
co^^  +  ycoG  -  co^ 


=  8oIm<! 


CO 


V^(l  +  Xo)-®^(l  +  Xoo)  +  7®cj(l  +  Xoo) 


2  2 

co^  -CO  -ycoG 


2  2 
COf  -CO 


-i2 


2  2 
+  CO  G 


:  8oCOG<i 


(l  +  Xcx))®f  “(l  +  Xoo)®  ~^t  (l  +  Xo)  +  ®  (l  +  Xoo) 


2  2 
COf  -CO 


2  2 
+  CO  G 


:  -8oCOG 


®^^(Xo-Xoo) 


2  2 
CO^  -CO 


2  2 
+  CO  G 


Then  the  loss  tangent  is  given  by 


tan  5  =  -- 


80COG 

Im8(co) 

ra^^lXo-Xoo) 

co^^  —  co^ 

']"  +  toV 

Re8(co) 

£0^ 

_®f^(l  +  Xo)-®^l 

(l  +  Xoo) 

V-co^]  +  coV(l  +  Xoo) 

1  1 
<N 

3 

1 

3“ 

2 

^  2  2 
+  CO  a 

=  COGi=^ 


®f^(Xo-Xoo) 


(1  +  Xo)-®  (1  +  Xoo)  C0/-C0"  +coV(l  +  x„) 
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The  expressions  above  are  valid  from  do  to  near-IR  frequenoies,  but  for  miorowave  applioations 
they  are  useful  primarily  at  low  frequencies.  It  is  therefore  convenient  to  take  the  following  limit: 

co^O^tanS^coa 

®r(l  +  Xo) 

The  problem  with  this  expression  is  that,  unlike  a  true  resistance,  this  quantity  depends  linearly 
on  frequency. 

In  order  to  make  contact  with  our  NLTL,  we  look  at  the  capacitance  of  a  BST  layer  with  area  A 
and  thickness  h  and  this  kind  of  dielectric  constant,  and  compute  the  corresponding  impedance  of 
the  layer: 

^  ,  .A  +  +  +  ^ 

C  =  8(co)-  =  8o  - 2 - ^ - 2 - h: 

n  cOf  +  ycoG-co  n 

^1  1  +  /COG  -  co^  h 

Z  =- -  =  ^ - < - - - - r - 1  — 

j(X)C  7CO80  ^  I- j  ^  ^  ^  _  (^2  j  ^  j  ^  ^  J  ^ 


A  low-frequency  expansion  of  Z  gives: 


Z  = 


1 

7C08o(l  +  Xo) 


A 


1  ^  I  f  /  “  r-  Yl  ^  + 

- 7 - \ J — - y  ^ - I 

7C08o(l  +  Xo)  ^  ^  co/  1  +  XojJ 


1  ^  L  I  y  to  _  1  o  Xq-X^  d 

7C08o(l  +  Xo)  CO, 2  I  1  +  Xo  Jj  7COSo(l  +  Xo)^  aoco,^  (l  +  Xo)^  ^ 

1 

= - hr 

J^Cq 
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and  thus  it  is  clear  that  Z  is  the  impedance  of  a  eapaeitance  Co  and  a  resistance  r  connected  in 
series.  This  justifies  our  cireuit  topology,  whieh  is  gratifying,  but  more  importantly  it  gives  us  a 
formula  for  r: 


tanS  =  COG- 


Xo-X. 


r  =  lim 


®i^(l  +  Xo) 
d  1 
C08o(l  +  Xo) 


CJ  Xo-Xoo  _  1  Xo-Xoo 

- 7 - T  - 1 - \  — 77 - 7 

(1  +  Xo )  (1  +  Xo )  to/  (1  +  Xo ) 

tanS  ^  Z  =  — — (l+  /tanS) 


tan  5 

= - 

coCq 


This  simple  result  almost  gets  us  where  we  want  to  be.  However,  in  order  for  r  to  be  constant  we 
need  the  linear  frequency  dependence  of  tan  5  alluded  to  above.  The  variation  of  tan  5  observed 
in  experiment  is  not  encouraging:  the  following  plots  (see  figure  8)  are  given  in  the  paper  by 
Vorobiev  and  Gevorgian  {10)\ 
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Frequency  (GHz) 


Figure  8.  Loss  tangent  of  a  BST  varaetor  vs.  frequeney  at  different 
temperatures  without  (a)  and  under  20  V  de  bias  (b).  The 
dotted  line  represents  the  tan  d  ~  power  law. 
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1/3 

These  authors  predict  a  peculiar  dependence  tan  co  for  the  loss  tangent,  which  is  not  very 
useful  for  our  applications.  Since  the  data  used  to  make  these  plots  was  unavailable,  the  crude 
approach  of  connecting  the  points  co  =  0,  tan  5  =  0,  and  co  =  40  GHz,  tan  5  =  0.07  with  a  straight 
line  was  used,  and  its  slope,  along  with  a  value  of  Co  appropriate  for  our  structure,  was  inserted 

in  the  formula  r  =  Using  the  endpoint  data  tanS  =  .07  at /=  40  GHz  gives 

COCg 


h 


r  = 


(.07) 


A  27r(4xl0^°s“')(8. 85x10“^^  f  cm“^)(5000) 
h 


1  f  =  1  Ohm  ' -s  ^  r  = 


(.07) 


A  27r(4xl0^°s“')(8.85xl0“^^  Ohm“^-s  cm“')(5000) 
h  1  (.07) 


A 


cm  27r(4xl0'‘^)(8.85xl0“^^cm“^)(5000) 


Ohm 


In  all  the  simulations  to  date.  Brown’s  numbers  have  been  used  for  the  line  dimensions.  In  his 

2 

notation,  A=wd=  (2.5  cm)(.3cm)  =  .75  cm  ,  A  =  .55  cm,  which  gives 


.55  cm  (-0^) 

.75  cm^  27r(4xI0'°)(8.85xI0“^^  cm“^)(5000) 


Ohm 


_ _ 

(6.2832)(.75)(4xl0'“)(8.85xl0“'‘‘)(5000) 


Ohm  =  .00046  Ohm 


As  better  data  is  obtained,  this  approach  can  be  refined.  For  now,  however,  the  value  obtained  is 
suitable  for  simulations. 


5.  Realistic  Antenna  Coupling  Via  an  RC  Circuit 


The  code  written  in  Mathematica  mimics  that  of  both  Brown  and  Ikezi,  i.e.,  it  is  based  on  a 
simple  series-parallel  circuit  model  for  the  NLTL  (see  figure  9). 
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For  the  applications  under  study,  this  transmission  line  is  excited  by  transient  signals,  i.e.,  pulses. 
Because  pulse  shapes  are  nonperiodic  functions,  there  is  a  potential  for  a  lossless  transmission 
line  to  support  signals  with  dc  levels.  However,  a  real  antenna  cannot  pass  a  dc  level,  so  that  this 
problem  is  clearly  mathematical  rather  than  physical.  To  represent  the  antenna  more  realistically, 
a  series  capacitance  can  be  included  in  the  load,  so  the  circuit  becomes  (see  figure  10). 


The  inductor  voltage  drops  are  related  to  the  loop  currents  as  stated  previously.  In  particular,  at 
the  last  stage  of  the  transmission  line  we  still  have 

dif 


L 


‘■N 


dt 


This  last  stage  sees  a  voltage  divider  described  by  the  equation 


T  _J/  D 

'^L  Hn 

dt 


=  C“V^_1  -(C“'  +n^_i  -{r  +  R^)i 
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We  can  use  these  expressions  to  eliminate  the  voltages  entirely,  which  gives  us  a  set  of  network 
equations  in  terms  of  and  alone: 


L^  =  Vs{t)-RJ^-T[q^-q2\-r[i^-i2) 


L =  r [q^_2  - q„.i ]  - r [^„_i  -qn\  +  r  (z„_2  - 2z„_i  +  ) 
^  ^  =  r  [q„_^  -q„]-T[q„-  q„^^  ]  +  r  -  2i„  +  ) 

L  =  r  -  g„+i  ]  -  r  ] + '^  (4  -  24+1 + 4+2 ) 


^  '  “^[^V-2  ^V-l]  ^Ar]  +  ^(iiV-2  2/jy_i+/^) 

^  ^  =  <2“ Vat-i  -  (C“'  +  +  a  (^^_1  -  )^  +  n^_i  -  (r  +  )  W 

Then  the  set  of  network  equations  is  completed  by 

.  _  dq„ 

"  dt 
.  _  dq^^ 
dt 

and 

It  is  no  longer  convenient  to  nondimensionalize  these  equations  because  there  are  two  different 
capacitors  in  the  system. 


6.  Description  of  Mathematica  Code 


The  circuit  equations  are  easily  converted  to  Mathematica™  form,  which  has  a  universal 
numerical  differential  equation  solver  called  NDSolve.  The  only  challenging  problem  about 
setting  up  these  equations  is  to  generate  the  code  for  an  arbitrary  number  of  layers.  This  can  be 
done  piecemeal  by  separating  the  equations  into  starting,  intermediate,  and  ending  segments. 
Each  of  these  segments  involves  a  set  of  equations,  as  shown  above: 
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Starting  segment: 


dq^  _  . 


dt 


=  h 


Intermediate  segment: 


L  =  r [q„.2  -  dn-i ]  - r [qn-i  -qn]  +  r{in-2  - \-i  +  i„ ) 


dqn-\  _  ■ 
dt 


n-\ 


di„ 


^ =  r [q„-i  -q„]-T[q„- q„^^ ]  +  r - 2i„  +  ) 


dq„ 

dt 


= 


L  =  ^[q„-  dn+l  ]  -  r  [qn+l  -  qn+2  2z„+i  +  4+2  ) 


dq 


n+\ 


dt 


=  i 


n+\ 


Ending  segment: 


dt 


C  +Ci  +n^_i  {r  +  R^^i 

dq^ 


N 


dt 


=  i 


N 


The  resulting  eode  looks  as  follows: 

eqs=Table [ {  cd[i] ' [tau] □ (l+LM  (s[i- 

1]  [tau]^2+s[i]  [tau] '^2+3  [i+1]  [tau] ''2-s  [i-1]  [tau]s[i]  [tau]- 
s [i] [tau] s [i+1] [tau] -s [i-1] [tau] s [i+1] [tau] ) ) (s [i-1] [tau] - 
2s [i] [tau] +s [i+1] [tau])+Qc  (cd[i-l] [tau]- 

2cd[i] [tau] +cd[i+l] [tau] ) , s [i] ' [tau] □cd[i] [tau] } , {i, 2 , NN-1 } ] ; 
eqst={cd[l] '  [tau] DUO [tau] - (s [1]  [tau] -s [2]  [tau] ) -LM  (s  [1]  [tau] - 
s[2]  [tau])''3-Qc  (cd[l]  [tau]-cd[2]  [tau]),s[l]'  [tau]ncd[l]  [tau]  }  ; 
eqfn={cd[NN] '  [tau] □ (s [NN-1] [tau] -s [NN] [tau] ) +LM  (s [NN-1] [tau] - 
s  [NN]  [tau] )  ^3-gaitima  s  [NN]  [tau] +Qc  cd  [NN-1]  [tau] - 
(Qc+QL) cd[NN] [tau] , s [NN] '  [tau] □cd[NN] [tau] } ; 
eqp=Prepend[eqs, eqst] ; 
equ=Flatten [Append [eqp, eqfn] ] ; 

init=Flatten [Table [{s[i] [0]n0, cd[i] [0] DO } , {i, 1, NN} ] ] ; 
NDSL=Flatten [Prepend [init, equ] ] ; 
vlist=Flatten [Table [ { s [i] , cd[i] } , {i, NN} ] ] ; 
sol=NDSolve [NDSL, vlist , { tau, 0 , 500 } , MaxSteps->Inf inity] ; 
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In  this  code,  the  function  Table  generates  an  arbitrary  number  of  intermediate  equations,  which 
describe  the  interior  of  the  transmission  line.  This  set  of  equations  (the  intermediate  segment)  is 
concatenated  with  the  starting  and  ending  segments,  and  the  full  set  is  solved  by  the  numerical 
routine.  The  result  is  a  list  of  currents  and  charges  for  the  cireuit  elements,  which  can  be 
examined  and  plotted. 

The  translation  of  this  mathematieal  routine  into  a  MatLab™  eode  will  be  straightforward, 
provided  that  MatLab™  has  utilities  which  can  generate  equation  lists  for  arbitrary  numbers  of 
layers.  Implementing  this  transfer  will  be  an  important  deliverable  when  work  on  this  project 
resumes. 


7.  Results  of  Preliminary  Simulations 


To  date,  the  numerical  code  has  been  exercised  in  ways  that  bring  out  the  physics  of  the 
frequency  translation  rather  than  parameter  extraetion;  this  focus  lends  eredenee  to  the 
correctness  of  the  code  results,  which  can  be  nonintuitive.  Indeed,  mueh  of  the  value  of 
numerical  work  lies  in  ereating  intuition  where  none  existed  before. 

Simulations  were  performed  using  a  series  of  pulses  based  on  Ikezi’s  normalized  trapezoidal 
pulse  shape.  I  chose  this  pulse  shape  beeause  it  ineorporates  simply  three  numbers  that 
characterize  any  pulse:  a  rise  time,  a  pulse  “width”,  and  a  fall  time.  Thus,  eonsider  a  pulse  of 
almost  triangular  shape  (see  figure  11). 


Figure  11.  Ikezi  normalized  short  test  pulse  shape;  time  scale 
in  NLTL  low-frequency  single-stage  periods. 


20 


For  a  lossless  NLTL,  the  time  seale  is  measured  by  the  period  of  the  LC  stages.  In  units  of  this 
period,  the  rise  time  of  this  pulse  is  2  periods,  its  duration  1  period,  and  its  fall  time  5  periods. 
Let  us  ehoose  L  and  C  sueh  that  the  period  was  13  ns,  and  pulse  the  line  with  a  50  V  signal, 
whose  total  pulse  length  is  104  ns.  This  low  voltage  should  cause  the  line  to  respond  like  a  linear 
filter.  Plotting  the  time  dependence  of  the  voltage  at  the  antenna  end  of  the  line  gives  (see  figure 
12). 


Figure  12.  Low  voltage  response  to  figure  1 1  test  pulse  at  end  of  NLTL. 


The  oscillatory  tails  following  the  pulse  front  exhibit  typical  dispersive  behavior  of  the  line. 
When  the  pulse  voltage  is  increased  to  6  kV,  the  following  plot  results  (see  figure  13),  which 
reveals  the  evolution  of  solitons. 


Figure  13.  Fligti  voltage  response  to  figure  1 1  test  pulse  at  end  of  NLTL;  note  solitons. 
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Increasing  the  50  V  pulse  width  to  10  periods  (see  figure  14). 


Figure  14.  Ikezi  normalized  long  test  pulse  shape. 
Gives  a  linear  response  of  the  form  (see  figure  15). 


Figure  15.  Low  voltage  response  to  figure  14  test  pulse  at  end  of  NLTL. 

The  longer  duration  has  separated  the  rise-time  regions,  produeing  spikes;  these  are  not  solitons, 
however,  since  the  voltage  is  too  low.  Increasing  the  voltage  to  5  kV  shows  multiple  soliton 
generation.  This  suggests  that  long-duration  pulses  are  better  than  short-duration  ones  for 
produeing  multiple  solitons  (see  figure  16). 
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Figure  16.  Fligh  voltage  response  to  figure  14  test  pulse  at  end  of  NLTL;  note  solitons. 


8.  Movies  of  NLTL  Interior 


There  are  a  number  of  features  of  these  plots  that  require  explanation.  To  this  end,  I  have  used  a 
very  handy  feature  of  the  Mathematica™  code  to  examine  the  interior  of  the  NLTL  as  a  pulse 
propagates  within  it.  Because  they  exhibit  only  the  output  voltage,  these  plots  treat  the  device  as 
a  “black  box”,  whose  mysterious  generation  of  alternating  “bursts”  of  pulses  is  hard  to 
understood.  However,  because  the  code  simultaneously  evaluates  voltages  and  currents  at  every 
node,  i.e.,  at  every  layer  of  the  structure,  these  nodal  quantities  can  be  exhibited  as  “maps”  of  the 
voltage/current  profile  at  each  time  step,  which  can  be  displayed  as  “movies”  of  the  line  interior. 

Multiple  movies  were  made  to  examine  the  overall  propagation  of  signals  in  the  line.  Two  of 
these  movies  are  included  in  this  report,  on  a  companion  CDs  for  the  hard  copies;  the  data  files 
they  are  based  on  are  very  large  (-.5  Gb).  Being  able  to  examine  the  interior  immediately  reveals 
the  origin  of  the  bursts:  the  initial  drive  pulse  propagates  to  the  end  of  the  line  with  a  certain 
average  velocity,  disintegrating  as  it  propagates.  Once  it  reaches  the  end  of  the  line,  it  is 
reflected.  In  the  simulations  I  have  made  to  date  the  antenna  impedance  is  high,  which  makes  the 
end  of  the  line  look  like  an  open  circuit.  This  causes  the  reflected  signal  to  add  to  the  incident 
signal,  doubling  the  total  amplitude.  The  first  pulse  burst  marks  the  arrival  of  the  distorted  pulse 
at  the  end.  The  pulse  then  returns  to  the  source,  where  the  voltage  is  fixed,  which  makes  the 
source  look  like  a  short  circuit;  this  causes  the  reflected  voltage  to  become  negative.  The 
negative  reflected  pulse  then  returns  to  the  far  end  of  the  line,  producing  the  negative-going 
second  burst,  which  is  once  again  doubled  in  magnitude.  This  cycle  repeats  until  the  pulse  has 
given  up  all  its  energy  to  the  load,  in  small  increments  because  the  load  impedance  is  high. 
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The  first  movie  is  of  a  long  (120  ns)  low-voltage  (50  V)  pulse,  whieh  ehanges  shape  as  it 
propagates  due  purely  to  dispersion.  Beeause  the  propagation  veloeity  of  a  sinusoidal  signal  in  a 
dispersive  medium  depends  on  frequeney,  the  high-frequeney  eomponents,  whieh  eontribute  to 
the  rapid  rise  and  fall  times  of  the  pulse,  “dephase”  during  the  propagation  and  produee  the 
oseillations  shown.  Note  that  the  frequeney  speetrum  tends  to  redistribute  energy  to  lower 
frequeneies  and  not  higher  ones.  The  seeond  movie  is  of  a  pulse  with  the  same  duration  but  mueh 
higher  voltage  (5  kV).  The  high-voltage  pulse  disintegrates  into  what  appear  to  be  solitons  before 
reaehing  the  end  of  the  line.  However,  their  shape  is  not  stable  -  there  is  an  unexpeeted 
oseillation  in  the  soliton  height.  Compare  these  pulses  with  the  superfieially  similar  “ripples”  in 
the  low-  voltage  pulse,  whieh  are  more  or  less  stationary  in  shape.  The  origin  of  the  high-voltage 
shape  instability  is  unelear;  while  it  is  tempting  to  identify  it  with  the  “modulation  instability”  of 
soliton  trains  discussed  in  the  literature,  it  could  also  be  related  to  the  discrete,  i.e.,  “circuit”, 
nature  of  the  model.  Although  this  effect  is  small,  and  hence  of  no  engineering  importance,  it  is 
scientifically  interesting,  and  cautionary  given  the  tendency  for  any  kind  of  pulse  compression  or 
amplification  to  be  referred  in  the  literature  to  as  “soliton  generation”. 


9.  Systems  Issues 


Although  this  modeling  effort  was  preliminary,  it  is  useful  to  bound  the  behavior  of  the  NLTL 
based  on  optimal  predictions  of  its  performance.  These  predictions  can  only  be  based  on  the 
information  available  now,  but  they  serve  as  a  guide  for  future  efforts. 

1 .  Maximum  output  power.  Clearly  we  can  use  the  breakdown  field  of  the  ferroelectric  to 
estimate  this  quantity.  In  BST  ceramics  this  field  can  approach  750  kV/cm.  Assume  that 
the  line  is  a  BST/alumina-filled  parallel-plate  waveguide,  which  propagates  a  TEM  wave 
with  electric  and  magnetic  fields  E  and  H  perpendicular  and  parallel  to  the  conducting 
walls  of  the  guide  (top  and  bottom),  respectively.  At  high  power  levels  the  dielectric 
constant  of  the  BST  saturates,  dropping  from  5500  for  a  very  good  sample  to  about  120. 
Then  for  a  uniform  line  made  of  lossless  material  the  characteristic  impedance  at  large 
signal  levels  is  =  377  /  Vs  =  377  /  Vl20  =  34.4  ohms .  Since  the  power  flux  is  given  by 
the  Poynting  vector,  the  power  flux  at  breakdown  is 

2 

p  =  fi’ExH  =  EH  =  Z~^E^  ={3AA  Q)"^ (7.5x10^  V/cm) 

=  16.3  GWW 

Since  the  lines  we  fabricated  had  cross-sections  of  2.5  x  .3  =  .75  cm^,  one  of  our  lines 
could  handle  '-'12  GW.  Note,  however,  that  the  inevitable  impedance  mismatch  with  the 
antenna  due  to  the  nonlinearity  will  reduce  this,  probably  by  about  50%,  so  -b  GW  is  a 
more  realistic  number. 
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2.  Maximum  frequency:  this  is  hard  to  specify  without  referring  to  geometric,  i.e., 
manufacturing  issues.  The  modeling  we  did  was  based  on  a  formula  Ikezi  quotes  in  one  of 
his  General  Atomics  papers  (Ikezi  et  ah,  Appl.  Phys.  Lett.  58,  986  (1991)),  derived  from 
simulation  data.  The  structure  he  starts  with  is  a  parallel-plate  waveguide  with  BST  sheets 
in  it,  i.e.,  the  one  described  above.  His  initial  design  numbers  involve  the  linear  modes  of 
this  structure,  which  is  in  essence  a  band-stop  fdter.  As  such,  it  allows  frequencies  to 
propagate  that  are  below  a  certain  maximum  frequency  -  he  calls  it  co^  -  and  blocks  signals 
with  higher  frequencies  within  a  “stop  band”.  His  nonlinear  simulations  reveal  that  this 
maximum  frequency  also  determines  the  oscillation  frequency  cOosc  from  the  empirical 
formula 

CO  ( 

^  =  1.25  +  1.51ogio  — 

J 

where  the  field  E  is  the  maximum  field  in  the  line,  and  Eq  is  a  material-dependent  constant 
with  the  dimensions  of  field;  for  BST  this  number  is  around  10  kV/cm.  This  frequency  co^ 
depends  mostly  on  the  widths  of  the  ferroelectric  and  dielectric  slabs,  and  their  linear 
dielectric  constants.  If  we  assume  the  ferroelectric  is  BST6-4,  i.e.,  with  60%  barium  and 
40%  strontium,  s  for  his  material  is  5500  at  room  temperature,  i.e.,  20  °C  (this  is  an 
optimistic  number);  at  lower  temperatures  it  is  higher  (the  Curie  temperature  is  0  °C),  but 
clearly  room  temperature  is  appropriate  for  our  application.  Given  the  volume  constraints 
imposed  by  the  latter,  the  optimal  transmission  line  length  is  probably  about  15  cm.  To 
increase  co^,  we  want  to  make  the  layers  as  thin  as  possible.  However,  we  found  that  layers 
thinner  than  .1  mm  are  too  fragile  to  handle,  so  let  us  assume  the  ferroelectric  layers  are  .1 
mm  thick.  If  we  have  15  cm  to  fit  the  transmission  line  into  the  weapon,  and  we  want  100 
layers  of  each  kind  of  material,  this  implies  a  total  thickness  of  ferroelectric  of  1  cm, 
leaving  14  cm  for  dielectric,  or  1.4  mm  per  layer.  Then  using  the  formula  above  along  with 
Ikezi’s  linear-design  formula  (equation  3  of  the  Applied  Physics  paper  mentioned  above), 
we  find  the  resulting  bandstop  fdter  corresponds  to  an  cOogc  of  3.4  GHz. 

3.  Pulse  repetition  frequency.  This  can  be  estimated  based  on  the  observation  that  in  this  type 

of  line,  impedance  mismatch  with  the  antenna  is  inevitable  due  to  the  nonlinearity.  This  in 

turn  gives  rise  to  reflections,  which  return  to  the  source  and  interfere  with  subsequent 

pulses.  The  safe  thing  to  do  is  to  wait  until  all  these  have  died  away.  If  we  define  “died 

away”  to  mean  that  only  1%  of  the  original  power  is  left  after  N  reflections,  and  the 

reflected  power  is  50%  =  .5  of  the  incident  power  for  each  reflection,  we  are  saying  that 

(.5)^=  .01,  or  A  =  between  6  and  7.  Now,  the  maximum  propagation  velocity  of  a  pulse  in 

the  transmission  line  is  the  velocity  of  light  for  the  reduced  (saturated)  high-field  dielectric 

constant,  which  is  120  for  BST  at  high  fields.  This  velocity  is  2.7  x  10  cm/s,  so  that  for  a 

—8 

15  cm  line  we  need  7  round  trips  or  a  distance  of  210  cm,  which  takes  7.7  x  10  s  =  77  ns. 
Then  the  PRF  is  the  reciprocal  of  this,  or  13  MHz. 
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Now,  suppose  we  have  a  train  of  30  pulses  coming  at  a  rate  of  3.4  GHz.  Then  the  time 
between  pulses  is  .29  ns  and  hence,  the  total  duration  of  the  pulse  train  is  30  x  .29  =  8.8  ns. 
This  implies  a  duty  cycle  of  1 1%.  Then  the  time  it  takes  for  the  first  pulse  to  hit  the  antenna 
is  5.6  ns,  and  the  time  to  return  to  the  source  is  11.1  ns.  Since  this  is  greater  than  8.8  ns,  all 
the  pulses  in  the  train  will  be  transmitted  to  the  antenna  before  the  first  reflected  pulse 
returns  to  the  source.  Since  the  power  is  6  GW,  over  a  period  of  8.8  ns  the  pulse  train 
delivers  53  J  to  be  radiated  by  the  antenna.  If  the  pulser  can  supply  a  kilojoule  of  energy, 

19  of  these  pulse  trains  can  be  generated  for  a  total  dwell  time  of  1.46  ps. 


Summary: 

Max  power: 

6  GW 

Max  frequency 

3.4  GHz 

PRF 

13  MHz 

Time  between  pulses 

290  ps 

No.  of  cycles  of  3.4  GHz/train 

30 

Duration  of  pulse  train 

8.8  ns 

Energy/pulse 

53  J 

No.  of  pulse  trains  in  dwell 

19 

Dwell  time 

1 .46  ps 

10.  Conclusions 


The  results  of  this  work  validate  the  previous  work  of  Ikezi  and,  to  a  lesser  extent.  Brown,  and 
show  that  we  now  have  a  useful  numerical  tool  for  studying  NLTL  and  guiding  efforts  to 
fabricate  them.  Future  work  will  center  on  parameter  studies  for  designing  and  optimizing  the 
performance  of  these  devices  as  HPM  sources. 
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Appendix.  Brown-Ikezi  Network  Equations 


The  circuit  is  defined  using  Kirchhoff  s  law  in  the  usual  way,  based  on  loop  currents.  Source  and 
load  resistances  are  introduced  to  simulate  coupling  to  the  PFN  and  antenna  (see  figure  A-1). 


The  inductor  voltage  drops  are  related  to  the  loop  currents  as  follows: 

L^=ys(t)-R,A-y, 


T^lnzL-y  _y 
^  -  r  „_2  y„_i 


L^  =  V„_,-V, 


L 


dt 

di 

dt 

di, 


n+\ 


dt 


=  v„-v„,, 


T  ^  =  Y  -Y 

dt  ^  ^ 


For  the  shunt  branches  we  introduce  the  nonlinear  function  T  [g]  =  C  Q  +  aQ  ,  where  C  is  the 
linear  term  in  the  capacitance  and  a  measures  the  circuit  nonlinearity.  As  discussed  in  the  text, 
we  introduce  “partial  charges”  q„  related  to  the  physical  currents  by  .  This  gives  the 


dt 


following  expressions  for  the  node  voltages: 
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Vo=Vs{t)-RJ, 

Vi  =  T[q^-q2\  +  r{i^-i2) 

Vn-I  =  r [qn-i  -qn\  +  - i„ ) 

K+l  =  r  [qn+l  -  qn+2  ]  +  r  {in+l  “  4+2  ) 

=  r  + r  ) 

VL=RLk 

We  can  use  these  expressions  to  eliminate  the  voltages  entirely,  which  gives  us  a  set  of  network 
equations  in  terms  of  and  z„  alone: 

=  ^5  {t)-Rinh  -r[^i -q2\-r{h -h) 


L 

L 

L 


=  r  [qn-2  -  qn-l  ]  -  r  [^„_1  -  ^„  ]  +  r{i„_2  -  2z„_i  +  Z„  ) 
^  =  r [q„_^  -q„]-Y[q„- q„^^ ]  +  r (z„_i  - 2i„  +  z„+i ) 

=  r  -  ^„+i  ]  -  r  ]  +  r  (z„  -  2z„+i  +  z„+2 ) 


^  ]  +  ''(%-!  2z^+4) 

^  ^  =  r  [^iv  -  ] + ''W  -  (^i  + '')  4 

It  is  useful  to  further  reduce  the  terms  involving  the  function  Y\Q\  as  follows: 

r  Un-X  -qnV^Un-  ^«+l  ]  =  {^«-l  “  1  +  {^«  “  ^«+l  1  “  “  ^«+l  1^ 

=  ( ^«-i  -  2^„  +  qn+l )  +  «  [{qn-l  -q„?-a{q„-  q„^i ] 

=  [C^'  +ai{q„_i -q„f  +{q„_i-q„}  [q^  -^„+i}  +  {^„  “2^„  +^„+i) 

=  [C^^  “2^„  +^„+i) 
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which  we  rewrite  as  follows: 

n  [qn-i ,  qn  ^  qn+i  ]  =  1  +  +  q„^  +  q„^^  -  q„_^q„  -  q„q„^^  -  q„+^q„_i ) 

2 

no  [qn-i ,  =  1  +  aC  (^„_1  -  q„ ) 

^  r [^„_1  -q^]-Y[q„- ]  =  C“'n [g„_i , q„ , q„^^ ]  * [q^_^  - Iq^  +  q^^y ) 
Then  the  final  set  of  network  equations  is 

.  _dqn 
"  dt 

and 

^ ^  =  ^5  (t)  - RJy  - C“'no  [^1 , ^2 ] ' (^1  - ^2 )  - ^{h  - h ) 

L =  C“'n [^„_2 , q„_y ,q„y  {q„_2  - 2q„_y  +q„)  +  - 2i^_y  +  /„ ) 

^  ^  =  <^“^n  [q^_y ,  q„ ,  q„^y  ]  -[q„_y  -  2q„  +  q^^y )  +  r  {i^_y  -  2i„  +  i„^y ) 

L =  C“^n [q„ , q„^y , q„^2 ] ' (^«  “ ^dn+i  +  ^«+2 )  +  ^ (i„  - 2z„+i  +  4+2 ) 

^ ^  =  <^“^n [^^_i ,qN,qL]  '{qj^_y  - 2^^  +  ^^ )  +  r{i^_y  - 2i^  +  ) 

L^^^Q[qN,qL\'{qN-qL)+ri^-{RL+r)iL 
where  the  state  variables  are  ^2’  ■"  >  (h’  4’  ■"  ’  'v’  h] 
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2  1 

Dimensionless  form:  let  cOq  = -  ,  =  rC  ,  =  R^^C  ,  and  =  i?^C .  Then 

LC 


_2  d\ 


®o“^  =  n  [g„-2,  ^„-l,  ]  •  [dn-l  -  2^„-l  +  ^„  )  +  Tc  (4-2  -  24-1  +  4  ) 

®o‘^  ^  =  n  [g„_i ,  qn ,  qn+l  ]  •  {qn-l  “  2g„  +  ^^+1 )  +  (4-1  “  24  +  4+1 ) 

®o“^  =  n  ,  qn+l,qn+2  ]'{dn-  ^dn+l  +  <in+2  )  +  (4  “  24+1  +  4+2  ) 


®0  Jl  -^[dN-l’^N’dlYidN-l  2g^  +  )  + ''^C  (^'lV-1 

®0  [dN’dlVidN  “^Z,)  +  '''C%  ~{'^L  +'''c)4 

Setting  T  =  coot ,  cooTc  =  (a^rC  =  Qc  ,  coqT^^  =  &oRi„C  =  Q.„  ,  coqT^  =  cooi?^C  =  gives 

^  =  ®0<214  (^)  -  no  [^1 ,  ^2  ]  ■  ®0  (^1  -  ^2  )  -  (Qin  +  Qc  )  h  +  Qch 

=  n  [qn-2 .  ^«-l .  •  ®0  {.(ln-2  “  2^„-i  +  ^„  )  +  Gc  (4-2  “  24-1  +  4  ) 

^  =  n  [g„_i ,  ]  •  ©0  (g„_i  -  2q^  +  )  +  Qc  (/„_1  -  2/„  +  4^1 ) 

=  n  [g„ ,  g„+i ,  g„+2  ]  ■  ®0  (^«  -  2^„+l  +  qn^2  )  +  Qc  (4  -  24+1  +  4+2  ) 

-  ^\q  N-\’d  N ’d  l\’  N-\  ~'^dN  +^i)  +  2c(^^-l  “2%  +4) 

“7^  =  ^0  [^Ar,^i]  ’^0  (^V  “^l)  +  2c^V  “(Gz,  +Gc)4 
ax 

Finally,  setting  Vg  (t)  =  F^u(t)  ,  q„  (t)  =  CF^q„  (t)  ,  4  (t)  =  cOoCF„ri„  (t)  gives 
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u  T1 

77  "  (0  -  no  [^1,^2  ]  •  (‘^i  -  ‘;2 )  +  a«®oCF^ri2  -  (a„  +  gc )  “oCF^rii 

W  'M 

=  n  [q„-2, ,<;„]■  ®oC7«  {^n-i  - 2q„-i  +  q„ )  +  Qc^qCV^  (ri„_2  - 2ri„_i  +  ri„ ) 

■  =  n  [q„-i, , q„+i  ]  ’  ®oCF„  (q„_i  -  2q„  +  )  +  Qc^^CV^  (ri„_i  -  2ri„  +  ) 

d  n 

■  =  n[q 

H»  ^«+l»  *5h+2  ]  i^n  ^in+1  <;«+2 )  +  Qc^qCV^  (t1„  -  2ri„+i  +  ri„^2 ) 

dw 

77  "  ^ (<;Ar-i  - 2q^  +  )  +  Qc^^CV^  (ri^_i  - 2ri„  +  ) 

^O^Kn  -  ^0  [‘^iV»‘^i]’®on7i  )  +  0C®on7^7V  ~(8l  8c) ’^O^Kn^L 


or 

^7TI 

^  =  U  (t)  -Ho  [^1,^2  ]  •  (<;i  -  <;2  )  +  Qin^l  -  [Qin  +  Qc)^\ 

dx 

dv\ 

=  n  [q„_2 ,  ,<;„]•  (<;„-2  -  2q„_i  +<;„)+ Gc  (n„-2  -  2ti«-i + ) 

d'w 

^  =  n  [q„_i ,  q„ ,  ]  *  (q„_i  -  2q„  +  [^n-\  “  2ri„  +  ri„+i ) 

=  n  [q„ , q„+i,  q„+2 ]  ’  “  2q„^i  +  q„^2 )  +  Qc  {^n  “  2ri„+i  +  ri„+2 ) 

d'w 

)  +  2c  (^a^-i  “2rijy  +ri2^) 
d  Ti 

[‘^7V»‘^l]’(‘’Af  “‘^i)  +  2c^Af  “(2i  +2c)^Z, 


where 

n  ?  *7,1+1  ]  ~  1  “I" "^*71,  *7«+l  ?«?«+!  ?«+!?«-! 

2 

no[^„-i,^„]  =  i+A(q„-i-qJ 

3  2 

and  A  =  aC  Vg  measures  the  system  nonlinearity. 
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NO.  OF 

COPIES  ORGANIZATION 

1  ADMNSTR 
ELEC  DEFNS  TECHL  INFO  CTR 
ATTN  DTIC  OCP 

8725  JOHN  J  KINGMAN  RD  STE  0944 
FT  BEL  VOIR  VA  22060-6218 

1  HC  DARPA 

ATTN  IXO  S  WELBY 
3701  N  FAIRFAX  DR 
ARLINGTON  VA  22203-1714 

1  CD  OFC  OF  THE  SECY  OF  DEFNS 
ATTN  ODDRE  (R&AT) 

THE  PENTAGON 
WASHINGTON  DC  20301-3080 

1  HC  US  ARMY  RSRCH  DEV  AND  ENGRG 
CMND 

ARMAMENT  RSRCH  DEV  &  ENGRG 
CTR  ARMAMENT  ENGRG  & 
TECHNLGY  CTR 

ATTN  AMSRD  AAR  AEF  T  J  MATTS 
BLDG  305 

ABERDEEN  PROVING  GROUND  MD 
21005-5001 

1  HC  PM  TIMS,  PROFILER  (MMS-P) 
AN/TMQ-52 
ATTN  B  GRIFFIES 
BUILDING  563 
FT  MONMOUTH  NJ  07703 

1  HC  US  ARMY  ARDEC 

ATTN  AMSRD  AAR  ALE  P 
C  HAINES 
BLDG  25 

PICATINNY  ARSENAL  NJ  07806-5000 

1  HC  US  ARMY  ARDEC 

ATTN  AMSRD  AAR  AEMC 
A  MARSTON 
BLDG  61 S 

PICATINNY  ARSENAL  NJ  07806-5000 

3  HCS  US  ARMY  ARDEC 

ATTN  AMSRD  AAR  AEML 
S  GILMAN 

ATTN  AMSRD  AAR  ALP  F 
O  NGUYEN 

ATTN  AMSRD  AAR  ALP  F  S  SHRI 
BLDG  65  S 

PICATINNY  ARSENAL  NJ  07806-5000 

1  HC  US  ARMY  INFO  SYS  ENGRG  CMND 
ATTN  AMSEL  IE  TD  A  RIVERA 
FT  HUACHUCA  AZ  85613-5300 


NO.  OF 

COPIES  ORGANIZATION 

I  HC  COMMANDER 

US  ARMY  RDECOM 
ATTN  AMSRD  AMR 
WC  MCCORKLE 
5400  FOWLER  RD 

REDSTONE  ARSENAL  AL  35898-5000 

I  HC  US  GOVERNMENT  PRINT  OFF 

DEPOSITORY  RECEIVING  SECTION 
ATTN  MAIL  STOP  IDAD  J  TATE 
732  NORTH  CAPITOL  ST  NW 
WASHINGTON  DC  20402 

I  HC  DEPT  OF  PHYSICS  AND  GEOLOGY 
UNIV  OF  TEXAS  -  PAN  AMERICAN 
ATTN  S  TIDROW 
1201  W  UNIVERSITY  DR 
EDINBURG  TX  78541 

I  HC  US  ARMY  RSRCH  LAB 

ATTN  RDRLCIMG  T  LANDFRIED 
BLDG  4600 

ABERDEEN  PROVING  GROUND  MD 
21005-5066 

25  HCS  US  ARMY  RSRCH  LAB 
ATTN  IMNEALCHRR 
MAIL  &  RECORDS  MGMT 
ATTN  RDRL  CIM  L  TECHL  LIB 
ATTN  RDRL  CIM  P  TECHL  PUB 
ATTN  RDRL  SE  J  PELLEGRINO 
ATTN  RDRL  SED  E  D  KATSIS 
ATTN  RDRL  SED  E  S  HENRIQUEZ 
ATTN  RDRL  SED  M  J  TATUM 
ATTN  RDRL  SED  P  K  JONES 
ATTN  RDRL  SEE  I  S  SVENS  SON 
ATTN  RDRL  SEE  O  W  M  GOLDING 
ATTN  RDRL  SLR  E  E  VIVEIROS 
ATTN  RDRL  SLR  E  F  CROWNE 
(5  COPIES) 

ATTN  RDRL  SLR  E  R  DEL  ROSARIO 
ATTN  RDRL  SLR  J  MAIT 
ATTN  RDRL  SLR  L  D  POTREPKA 
ATTN  RDRL  SLR  L  J  PULSKAMP 
ATTN  RDRL  SLR  L  M  DUBEY 
ATTN  RDRL  SLR  L  R  POLCAWICH 
ATTN  RDRL  SLR  P  AMIRTHARAJ 
ATTN  RDRL  SLR  U  C  FAZI 
ATTN  RDRL  SES  M  PATTERSON 
ADELPHI  MD  20783-1197 

TOTAL:  40  (I  ELEC,  I  CD,  38  HCS) 
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